typeset using JPSJ.sty <ver.l.Ob> 



OS 



Superconductivity of Quasi-Two-Dimensional Tight-Binding Electrons in a 



^ ' Strong Magnetic Field 

^ I 

Ch I Mitake Miyazaki, Keita KiSHiGi and Yasumasa Hasegawa 

^ ! 

^ ; 

^ _ ^ Faculty of Science, Himeji Institute of Technology, 

' Kamigouri-cho, Akou-gun, Hyogo 678-1297, Japan 

, (Received ) 

o ■ 

V ■ 

^ ' We have investigated the transition temperature Tc{H) of superconductivity in quasi- two- 

^ ' dimensional (Q2D) tight-binding electrons in a strong magnetic field. When the magnetic field 

■ is parallel to 2D conducting plane, Tc{H) of the Q2D superconductor is shown to increase in an 

. oscillatory manner as the magnetic field becomes large and to reach Tc(0) in a strong magnetic 

^ . field limit for the spin-triplet superconductor. We consider the cases of on-site and nearest 

'"^ , sites attractive interaction, and calculate the magnetic field dependences of the transition tem- 

^ I perature for various types of symmetry. The first order transition from py-wave to p^-wave is 

^ ' shown to occur at ~ 35T when the magnetic field is parallel to the y direction, which will 
be observed in a triplet superconductor, Sr2Ru04. 

> ■ 
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, §1. introduction 

' Recently, it has been shown that the quasi-two dimensional (Q2D) superconductivity evolves from 

^ I the Ginzburg-Landau-Abrikosov-Gor 'kovB'i (GLAG) and Lawrence-DoniachilB& (LD) region to 

. the reentrant phase as a magnetic field becomes large,i'0^ which is derived by considering the 

O ■ quantum effect of electron motion in the presence of a strong magnetic field apphed parallel to 2D 



X 
J3 



conducting plane. The Q2D systems with layered structure, for example high Tc oxides, /3-(BEDT- 
TTF)2l3& and Sr2Ru04,i^ have a two-dimensional cylindrical Fermi surface with weak warping 
along the kz direction. When the magnetic field H is parallel to the conducting plane, the electrons 
move in orbits which are the intersection of the Fermi surface and a plane normal to in a semi- 
classical picture. Some orbits are open and the others are closed. The former gives a similar effect 
as QID superconductivity in a strong magnetic field. By neglecting the contribution 
from the closed orbits, Lebed and Yamajii^ calculated the mean field transition temperature of a 
Q2D superconductor in the parabolic band model and demonstrated the possibility of reentrant 
superconductivity. In a previous paper,& we calculated the mean field transition temperature 
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numerically by taking account of the eigenstates of three-dimensional tight-binding electrons in 
a strong magnetic field. In this formulation, we can treat QID and Q2D in the same manner by 
changing t{j/ta, where ta and ti, are hopping matrix elements in the conducting plane. In the previous 
calculation,B however, we considered only on-site attractive interaction for the Q2D systems. 

The high Tc oxides are thought to be Q2D d-wave superconductors and k-(BEDT- 
TTF)2Cu[N(CN)2]Br§'@lll'lli) is shown to be anisotropic Q2D with hne nodes of the gap. Rice 
and SigristEl) proposed that Sr2Ru04 could be a spin-triplet superconductor by analogy with ^He 
and experimental evidences for the triplet superconductivity have been observed. 

TheQ2D 

anisotropic superconductivity in a strong magnetic field is of great interest. 

In this paper, we assume the nearest site attractive interaction in tight-binding model, and we 
examine the transition temperature of anisotropic superconducting states in a strong magnetic 
field. 

§2. Model and Green Function of Tight-Binding Electrons in a Magnetic Field 

The tight-binding electrons in a magnetic field are described by the Hamiltonian (we take h, 
and the velocity of light to be 1), 

Ti. = Ti-o + Tiu 

no = -ta E e''-4<xCi,<x-i6 E 

— tc E ^ '^irr'^i.g- ~ ^ E ^i.a'^i,<^ 

(i,j)c,cr «,o" 

- E <y^^BHc\^ci^a (1) 

'Hu= E f^ii4,aCi,acJa'Ci,a'> (2) 

<i,j>,a,u' 

where and q^o- are creation and annihilation operators at the i site, jx is the chemical potential, 
a/iBH is the Zeeman energy for | (J,) spin (cr = -|-(— )), Uij is the interaction between electrons at 
i and j sites and 

eij = '^rAdi. (3) 

In the above, A is the vector potential and 0o is the flux quantum. We consider the Q2D system 
with the hopping matrix elements ta > tf, ^ t^. The interlayer hopping matrix element t^. is taken 
to be larger than transition temperature for H = 0, so that the fluctuation due to the low 
dimension is small and we take the mean field approximation in this paper. 

In the interaction Hamiltonian, eq.(2), we take the on-site interaction and nearest-site interaction 
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along each axis as 



U, 



Uo if Ti = Tj 

Us if Ti = Vj ±S 
otherwise 



(4) 



where 8 = x, y and z are unit vectors along a axis, b axis and c axis, respectively. The Fourier- 
transform of the interaction is obtained as 



U{q) = Uq + 2Ux cos{aqx) + 2Uy cos{bqy) 
+2Uz cos{cqx). 



(5) 



In this paper, the magnetic field H is parallel to the 6 axis. We use the vector potential A = 
{0,0,— Hx). Since the magnetic field considered in this paper is not extremely strong {(t)/(j)o 1, 
where = Hac is flux per unit cell) , the change of the electron number as a function of H may be 
very small in QID and Q2D case. This can be understood by noting that there are no closed orbits 
in the QID case and only small part of the orbits are closed in the Q2D case in the semi-classical 
picture. Therefore, the chemical potential /x is fixed for given ta, % and tc to give a quarter-filled 
band H = instead of electron number to be fixed. 

The noninteracting Hamiltonian is written as 



M_i 




V 

Mo y 
y* Ml 



V 



V 



(6) 



where 



Mji = —2taCOs[a{kx + nG)] — 2tbCOs{bky) 
—gubH — fi. 



'ickz 



G={G, 0,0) 



-,0,0 



(7) 

(8) 
(9) 

(10) 



Cl = {■■■, ct(fc - G), 4 (fc), 4 (fc + G), ■ ■ ■), 

'27r 

a • 

The size of the matrix in eq.(6) is g x g if (j)/(f)Q = p/q, where p and q arc mutually prime integers 
and infinite if (f>/(f>o is irrational. The summation in k should be done in a magnetic Brillouin zone, 
\kx\ < t^/{qo), \ky\ < n/b, and \kz\ < tt/c. The eigenvalues of the matrix in eq.(6) for fixed ky 
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have very rich structure as a function of the magnetic field if ~ ta as shown by Hofstadter. 
In the case of tc <C ta, which is the case for QID and Q2D systems, however, the energy is almost 
continuous except for the large gaps near the bottom and top of the energy. The gaps near the 
bottom and the top can be interpreted as the result of the Landau quantization in the closed orbit 
and the negligible gaps can be understood as the small probability of the magnetic breakdown in 
the open orbit. Since we consider the instability to the superconductivity, only the states near 
the Fermi surface with the energy |e — < T are important. For the electrons of less-than-half 
filled band, or the quarter-filled electrons, we can consider only the part of the matrix in eq.(6) 
which has, for example, \a{kx + nG)\ < (3/4)(7r/a). By this approximation (j)/4>o is not restricted 
to be a rational number. The Brillouin zone should be taken as \kx\ < G/2 instead of the magnetic 
Brillouin zone \kx\ < Tr/{qa) = G/{2p) for the rational magnetic field case {4>/(f)o = p/q)- The 
energy does not depend on in this approximation. This can be seen by changing ct{k + mG) to 
g-imcfc^gt^^ + mG), by which only {l,q) and {q, 1) elements in eq.(6) depend on kz but they are 
irrelevant for the energy near the Fermi level in the case of tc <^ ta- The effects of the closed orbits 
and the magnetic breakdown are taken into account correctly in this approximation as shown in 
Fig. 1. 
We get the energy 

^n,k,a = ^("•' ~ '^^b cos{bky) - a - n (11) 
and the eigenstates \^ain,k)) by numerically diagonalizing the matrix of size of the order of 
[(3/4).27r/G]x[(3/4).27r/G]. 



The creation operators ct{k + mG) is expressed by the creation operators of the eigenstates 
{^l{n,k)) as 

ct (fc + mG) = e'^^'^' <t>lM. n)^>l{n, k), (12) 

n 

where the coefficient (/)fc^(m,n) can be calculated numerically. Notice that (/)fc^(m,n) does not 
depend on ky and kz . 
The real space Green's function is given by 

G^{n,r'j,ioJi) = -j\re'^^- (T,Cr-„.(T)cj;, _^(0)^ 
^ ^ ^ (PkArn,n)(j)l^{m',n) 

k,nm,m' ^'^l ~ ^n,k,a 

^ ^i{r'^^ri)-k+i{m'r'^^mri)-G 

^ gi(m'-m)cfc,^ (13) 

where uji = (21 + l)7rT is the Matsubara frequency. 
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Fig. 1. Fermi surface in kx-ky plane and energy as a function of kx in Q2D. Three curves in the upper figure are 

the Fermi surface in the pianos of kz = tt/c, 7r/(2c) and 0, respectively. The lower figure shows the energy as a 
function of kx for four values of ky shown in the upper figure, representing no-orbit (A), closed orbit (B) and two 
open orbits (C) and (D) from the top. We neglect the region of \akx/Tr\ > 3/4. In the presence of magnetic field 
gaps open at kx = mG/2 with integer m =^ 0. 



The linearized gap equation for isotropic and anisotropic pairing in coordinate representation is 
obtained in the mean field approximation as 

A^^>{ri,rj) = UijY,^K^^>{ri,rj;r'^,r't)A^^,{r'^,r't), (14) 

s t 

where 

K^a'{ri,rj;r'g,r[) = Tj2Ga{ri,r'g,iu;i)Ga'{rj,r[,-iui). (15) 

§3. Spin-triplet superconductivity 

In this section, we consider the spin-triplet superconductivity. We consider only the equal-spin- 
pairing state in the following. The Zeeman effect does not play any important role in a strong 
magnetic field, as long as we neglect a change of the density of states of the up and down spin. 

For the spin-triplet case, we define "p^-wave" the order parameter as 

A^^(r)^ i [A^^{r,r + S)-A^^{r + S,r)] . (16) 

Here, the Fourier transformation of the order parameters is defined by 

A^^ll(q) = J (irA^|(r)e-'(«+^^)-(^+^/2)^ ^^7^ 
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and 

q,N 

where Qx is taken as —G/2 < Qx ^ G/2 and N is integer. The p^-wave state can be a solution of 
eq. (14) if Us < 0. When the system has tetragonal symmetry, px and py have the same coupling 
constant at H = and px + ipy will be stabilized at T < Tc, resulting in the nodeless gap in the 
Fermi surface. The magnetic field along the b axis lifts the degeneracy of px and py even in a week 
magnetic field as discussed by Agterberg® in the GL theory. 

The linearized gap equation, eq.(14), is written as a matrix equation for A^ ||(q). The transition 
temperature is obtained by the condition that the nontrivial solution exists. It can be shown that 
the maximum transition temperature is obtained for q = 0. The matrix equation is reduced to the 
equations for even N and odd N. For even N the linearized gap equation is obtained as 

1 -1 



k^ky I'nn' 



^^l - £n,k^,ky^ ^^l + ^n'-ko,-ky-] 

X E E '^™i-"^2,m3-m4'5mi+m3,-Af<^^^_^g iV_ iV^ 
mim2 m3m4 ^ ^ 

X 4>k^ ("T-i , n)(l)k^ (m2, n)(p^k^ (ms, n')(l)^k^ (m4, n') 
N N'' 

771-4 — "13 



COS 



— COS 



— + — ]G-S 

2 2) 

+ [mi + m3 + — + — \G-5 



(19) 



where 



n 



ll,2v{<lx) = E EE(^fcm)^ 

ka:,kyn,n' ni 

Xfpka^im - l,n)(j)k^{m - l',n) 
X(/>_fc^(-m - l,n')<t)-kS-fn - l',n') 

^ ~ f {^n,kx,ky^t) ~ f {^n' -kx-qx,-ky^) 



2{£n,kx,ky1 + 



-kyV 



(20) 



where f{e) is the Fermi distribution function and 7^^^ has the following forms for each order 
parameter: 



For odd N, we get 



sm[a{kx — mG)] T=px-wave 
s'm(bky) T=py-wa,ve 
sm{ckz) T=pz-wave 



(21) 



(22) 
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where 

kx ,ky n,n' 

X(pkxim-l, n)(j)kx - I', n) 

X'P-kxi-m - I - l,n)(p^kx{-m -I' - l,n') 

^ 1 ~ f {^n,kx,ky) ~ f i^n' ,-kx-qx,-ky) ^^g-j 
'^{^n,kx,ky ~^ ^n' ,—kx-qx,—ky) 

The transition hne is given by 1 — gUs = 0, where g is the maximum eigenvalue of the matrix 
n for even and odd A^. In this paper, we calculate the field dependence of g at low temperature 
instead of calculating the transition temperature. 

In the following, we consider the Q2D superconductor of the quarter-filled electrons with param- 
eters = 1) tc/ta = 0.05 and T/ta = 0.001. In Fig. 2, we plot the effective coupling constant 
g/go for each states as a function of aG/{2-K) = (j)/(pQ, where go is the effective coupling constant 
for tc = 0, which corresponds to that in the absence of magnetic field. For example, the magnetic 
field of 20 Tesla corresponds to 0/(/)o ~ 0.0036 for a = 3.871 and c = 12.71 (Sr2Ru04). It is 
found that g/go obtained by diagonalizing the even part of the matrix 11 increases in an oscillatory 
manner as a magnetic field becomes larger and reaches that in the absence of magnetic field, while 
g/go for the odd part becomes zero in the strong field limit. In the semi-classical treatment of 
the magnetic field, superconductivity is destroyed in a strong magnetic field, since the energies of 
electrons with wave numbers k and —k is not equal. However, as is seen in eq. (19), the instability 
to the superconductivity is caused by the pairs between the states (n, k) and (n', —k) with the 
same energy en,kx,ky ~ ^n' -kx-ky- On the other hands, the coefficient (/>fc^(m,re) becomes small in 
a weak magnetic field region but a lot of states within the energy range of T contribute to forming 
Cooper pairs, which will reproduce the GLAG result. 



We obtain g/go is larger in the p^^-wave state than in the py-wa,ve state for (p/(j)o ^ 0.0065 as shown 
in Fig. 2. The difference oi g/go can be understood as follows. Since the electron orbit at /cy ~ is 
open, we can use the linearized dispersion and we can treat it as a quasi-one dimensional system.!^ 
Then we obtain (/)fc^(m, n) ~ Jm-n i'2tc/vY{ky)G), where J is the Bessel function and vp{ky) is the 
/cy-dependent kx component of the Fermi velocity. & Therefore, for the larger vp{ky), the larger 
effect of the magnetic field is expected. On the other hand, the electron orbit at ky ~ 7r/26 in the 
quarter filled band is closed or has a small vp{ky) even if the orbit is open (see Fig. 1). In the 
Pa;-wave state the amplitude of the order parameter is largest at kx = /2a, which corresponds to 
/cy ~ in the Fermi surface of the quarter filled band. The order parameter in py-wave state is zero 
at ky = 0. Thus the quantum effect, which causes the increase of g/go in a strong magnetic field, 
is larger in the p^^-wave state than that in the py-wave state. This is a simple interpretation using 
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Fig. 2. (a) Effective coupling constant as a function of in the case of tf,/ta = 1.0, tc/ta ~ 0.05, T/ta = 0.001 

and H \\ y. Two lines for each state are obtained from the even and odd parts. The larger g/go are shown by thick 
lines, (b) The low field region. 



the Fermi surface in H = 0. In the case of H ^ 0, we can interpret in the same way as follows. By 
using the approximation of the linearized dispersion, we can write eq. (20) for the px and py wave 
as, 

sin^ a{kx — kF{ky)) 



^21,21' — X! 

1 - f{VF{ky){kx - k^{ky))) - f{VY{ky){kx + iVG " k^{ky))) 



s\v? bky 



V^{ky){2kx + NG - 2kY{ky)) 

deJN-2i [ sine] Jn-2V [ f!\r^ ^ ) ' (^4) 

\vFiky)G J \vF{ky)G J 

where k-p{ky) is the /c^-dependent Fermi wave number. The large contribution comes from the region 
of large vp{ky), i.e. ky ^ 0, where sin^ a{kx — k^^ky)) > sin^ bky. As a result, the g/go for ^^.-wave 
is larger than that for py-wave. In the tetragonal system with ta = tc and Ux = Uy ^ Uz, the 
Px state is realized at the transition line in the strong magnetic field along the b axis {y direction), 
and the py component will become finite at the second transition at lower temperature. 

In the weak magnetic field, we can apply the GL theory. As shown by AgterbergJ^l) only the 
Py-wave orders at Hc2 (or at the transition temperature Tc{H)) when H is along the b axis, and 
at the second transition the p^.-wave becomes finite. The difference of the transitions for the px- 
and Py- waves comes from the diff^erence of the GL parameters ki and K2 (ki > ^2), which are the 
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coefficients of the gradient terms (|D^A^|(r)p + |Dj,A^^(r)p) and {\DyAP''^{r)\'^ + |D^A^|(r)|2), 
respectively, where = Vj — eAj. 

At a first glance, the GL theory and our result obtained in the strong field may be considered 
to be incompatible, but it is reasonable. As discussed above the quantum effect of the magnetic 
field to increase g/go in sl strong magnetic field is stronger in the Px-w&ve than in the py-wave for 
H \\ b. In the same way the orbital frustration effect to reduce the transition temperature in a 
weak magnetic field is stronger in the Px-w&ve than in the Pj,-wave, resulting in the lower transition 
temperature in the Px-w&ve in a weak magnetic field. In the crossover region of the magnetic field 
relative value of g/go for the Px-wave and the py-wave will change. In Fig. 2(b), we see these 
features at (t)/^o ~ 0.0065 {H ~ 35T). 

§4. Spin-singlet superconductivity 

In this section, we consider the spin-singlet superconductivity. In the presence of the nearest-site 
interaction, A^i{r,r + S) = — Aji-(r, r -|- 5) = A'^^{r + 6,r) = — Aj^|(r-|-^, r). Therefore, the order 
parameters of "d^-wave" is defined by 

Afj(r) ^ ^[A^^ir,r + S) + A^^{r + d,r)] (25) 

The dx-wave or dy-wavc is realized only ifUx < and Uq = or Uy < and Uq = 0, respectively. In 
contrast to the spin-triplet superconductivity, if Uq 7^ 0, dx and dy are not self-consistent solutions 
of eq. (14) but the linear combination of A| J (r) , A| J (r ) and A|| = A| j (r, r) can be a self-consistent 
solution. If the system has a tetragonal symmetry (ta = h and Ux = Uy < 0), the d-wave (dx — dy) 
is realized at H = 0. However, dx — dy is not the self-consistent solution in the presence of the 
magnetic field along the b axis. The dz is realized as long as Uz < and the higher terms in tc/ta 
are neglected. In this section we consider the case that one of the interactions (Uq or Us) is negative 
and the others are zero, so that s or d^-wave is realized, for simplicity. Here, we write 

Ajj(r) = J2 e'(^+^^)-(^^+^/') A^^|;(q). (26) 
q,N 

By taking S = 0, we get the s-wave order parameter. 

The linearized gap equation is written as a matrix equation for A^|j^(q) as in the spin-triplet 
case. The matrix equation is separated into even N and odd N parts. We get 

where coupling constant A is C/q or Us and 

n|,.= EEE(7L)' 

kx,kyn,n' m- 

X(pk:,(m - l,n)(f)k^ (m - l',n) 
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X (t)_kx {-fn - I, n')(f)-kx {-m - I', n') 

^{^n,kx,ky,'\ ^n',-kx,-ky,l) 

where has the following forms for each order parameter: 



1 <S=s-wave 

cos[a(A;x- — m,G)] S=dx-weLve 
cos{bky) 
cos{ckz) 



S=dy-wa,ve 
«S=(iz-wave 



For odd N, we get 



where 



^ X] ^2l+l,2l'+lU^2l'+lU^ 



kx,ky n,n' 

yi(t)-kx{-m - I - l,n')(p-kx{-m - l' - l,n') 

^ 1 - f{^n,kx,ky^) - f{£n>-kx-ky,i) 



(28) 



(29) 



(30) 



(31) 



2(en,fc^,fcj,,t ~^ ^n' ,-kx,-ky,l) 

In the following, we neglect the Zeeman energy for simplicity to show the different behavior 
between each state obviously. We first calculate the effective coupling constant of s-wave super- 
conductivity in a strong magnetic field and examine its tc/ta dependence. 



In Fig. 3, wc plot the effective coupling constant as a function of (t)/(j)Q- The effective coupling 
constant reaches that for = as a magnetic field increases. We find that the effective coupling 
constant depends strongly on the hopping matrix elements between layers tc- As tc/ta becomes 
small, the oscillation of g/go becomes small, and the value of g/gQ increases in whole. Thus, 
reentrant behavior will be observed in weaker magnetic field in the superconductor with smaller 

We now study the effective coupling constant of Q2D anisotropic superconductor with tc/ta = 
0.05. In Fig. 4, we plot g/go obtained by each pairing state as a function of ^/(/>o- 



For each order parameter, g/go of the even part increases and that of the odd part decreases in 
a strong magnetic field limit. The effective coupling constant g/go in the dy-wave state is larger 
than that in the dx-wavc state. This can be understood as in the spin-triplet case. The dy-wave 
have a largest order parameter at ky = and its order parameter is zero at ky = 7r/26, while the 
order parameter of the dx-wave is largest at kx = {ky « 7r/26 for the quarter filled band) and 
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0.40 I ' ' ' ' ' ' ' 1 

0.000 0.010 0.020 0.030 0.040 

Fig. 3. Effective coupling constant of the s-wave pairing as a function of 0/0o in the case oltb/ta = 1-0, T/ta = 0.001 
and H \\ y. 



zero at kx = 7r/2a {ky f« 0). As a result, the recovery of g/go in the strong magnetic field due to 
the quantum effect in the open orbit occurs in lower magnetic field for the dy state than for the dx 
state as seen in Fig. 4(b). 

§5. Conclusion 

In this paper we study the anisotropic superconductivity of tight-binding electrons with hopping 
matrix elements ta = ^6 ^ tc as a function of the magnetic field parallel to the b axis. We calculate 
the energies and the eigenstates in the magnetic field by numerically diagonalizing the matrix and 
the effective coupling constant is calculated by using these values. The effects of both open and 
closed orbits are taken into account in our calculation. We consider the attractive interaction 
between electrons in the nearest sites along each axis to realize the anisotropic superconductivity. 
With these interactions singlet {dx, dy and dz) and triplet {px, Py and Pz) superconductivities 
are possible. The effective coupling constant, g/go, for these anisotropic superconductivities are 
calculated. It is obtained that g/go approaches to 1 as becomes large and it depends on the 
symmetry of the order parameter. As shown in Figs. 2 and 4, we find that the Px-wave {dy-wscve) 
state gives higher transition temperature than that of py-wave (dx-wave) in the strong magnetic 
field region, and both of them reach the transition temperature of zero magnetic field in a strong 
field limit. The first order transition from py-wavc to p-^-wave is predicted at H ^ 35T, which will 
be observed in the p-wave superconductor, Sr2Ru02. 
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Fig. 4. (a) Effective coupling constant as a function of (ji/cjio in the case of tb/ta = 1-0, tc/ta = 0.05, T/ta = 0.001 
and H \\ y. Two lines for each state are obtained from the even and odd parts. The larger g/go are shown by thick 
lines, (b) The low field region. 
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